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Abstract 

^ In this paper a new multiscale modeling technique is proposed. It 

relies on a recently introduced measure-theoretic approach, which allows 
to manage the microscopic and the macroscopic scale under a unique 
framework. In the resulting coupled model the two scales coexist and 
share information. This allows to perform numerical simulations in which 
the trajectories and the density of the particles affect each other. Crowd 
dynamics is the motivating application throughout the paper. 
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1 Introduction 

Modeling group dynamics of living systems, such as groups of animals or hu- 
man crowds, is a difficult task because one can only partially rely on the well- 

i 1 established theories of classical mechanics. Mathematical models must take into 

account several features of living matter: for example, individuals are not pas- 
sively dragged by external forces, instead they have a decision-based dynamics; 
they experience nonlocal interactions, since they are able to see even far group 
mates and make decisions consequently; interactions can be metric (i.e., with 

V.Q group mates less than a threshold apart) or topological (i.e., with a fixed num- 

ber of group mates no matter how far they are) pQ; interactions are strongly 
anisotropic because the subjects have a limited visual field, and mechanisms 
for collision avoidance are expected to be mainly directed toward group mates 
in front 9J; individuals are different from each other, each of them having for 
instance her/his own goal, reaction time, and maximal velocity. 
^ One of the most interesting consequences of these characteristics is the emer- 

J^j gence of self-organization. Individuals can deploy themselves to give rise to ap- 

parently ordered and coordinated configurations or patterns |21j . We cite, for 
example, clusters by starlings [TJ, lines by elephants, penguins, and lobsters, 
V-like formations by geese, lanes by pedestrians [HI [TS] . Actually such group 
configurations are not the result of a common decision made by the individuals 
or by a leader. Instead, they stem from simple rules followed by each individual, 
which takes into account the position/velocity of a few group mates. It is then 
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possible that a single individual does not even perceive the global structure of 
the group it is part of. 

In this paper we focus on crowd dynamics, which in the last few years has 
been the object of many mathematical models. In the microscopic (i.e., agent- 
based) approach pedestrians are considered individually. Models usually consist 
of a (large) system of ordinary differential equations, each of which describes the 
behavior of a single pedestrian [TU [T51 US [351 [35] . In the macroscopic approach 
pedestrians are instead described by means of their average density, which in 
most models obeys conservation or balance laws [21 El EH [19] . 

It is not fair to state that either approach is better for whatever problem. 
Rather, it is clear that a microscopic approach is advantageous when one wants 
to model differences among the individuals, random disturbance, or small envi- 
ronments. Moreover, it is the only reliable approach when one wants to track ex- 
actly the position of a few walkers. On the other hand, it may not be convenient 
to use a microscopic approach to model pedestrian flow in large environments, 
due to the high computational effort required. A macroscopic approach may 
be preferable to address optimization problems and analytical issues, as well as 
to handle experimental data. Nonetheless, despite self-organization phenomena 
are often visible only in large crowds [15] . they are a consequence of strategical 
behaviors developed by individual pedestrians. 

In [TU1 [5S1 ISO] we have extensively analyzed a measure-based modeling frame- 
work able to describe group behavior at both the microscopic and the macro- 
scopic scale. The key point is the reinterpretation of the classical conservation 
laws in terms of abstract mass measures, which are then specialized to singular 
Dirac measures for microscopic models and to absolutely continuous measures 
(w.r.t. Lebesgue, i.e., the volume measure) for macroscopic models. We have 
shown jTU] that the two scales may reproduce the same features of the group be- 
havior, thus providing a perfect matching between the results of the simulations 
for the microscopic and the macroscopic model in some test cases. This moti- 
vated the multiscale approach that we propose here. Such an approach allows 
to keep a macroscopic view without losing the right amount of "granularity" , 
which is crucial for the emergence of some self-organized patterns. Further- 
more, the proposed method allows to introduce in a macroscopic (averaged) 
context some microscopic effects, such as random disturbances or differences 
among the individuals, in a fully justifiable manner from both the physical and 
the mathematical perspective. In the model that we propose, microscopic and 
macroscopic scales coexist and continuously share information on the overall dy- 
namics. More precisely, the microscopic and the macroscopic part of the model 
track the trajectories of single pedestrians and the density of pedestrians, re- 
spectively, using the same evolution equation duly interpreted in the sense of 
measures. In this respect, the two scales are indivisible. This makes the differ- 
ence from other ways of understanding multiscale approaches in the literature. 
For example, in |31j a multiscale geometric technique is used to represent the cir- 
culatory system: one specific part of the network is accurately modeled in three 
dimensions, whereas the rest is described by means of lumped zero-dimensional 
models. This enables one to account for the whole circulatory network while 
keeping the complexity of the model under control. Multiscale methods can be 
implemented also at a numerical level in connection with domain decomposition 
(see e.g., [HI [35]), m order to compute the solution to a certain equation with 
different local accuracy. The general idea is to couple accurate but expensive 
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calculations, performed by a microscopic (e.g., particle-based) solver in small 
and inhomogeneous regions, with less accurate but also less expensive ones, 
performed by a macroscopic (e.g., continuum) solver in large and homogeneous 
regions. The two solvers usually exchange information at the interface of the 
respective regions. Another possibility (see e.g., [2D]) is to alternate the two 
scales for computing on the same system. In the resulting iterative algorithm, 
the output of the microscopic simulation is used as input for the macroscopic 
simulation and vice versa. A further way to understand the multiscale approach 
is through upscaling procedures. In this case, the ultimate goal is to pass from 
a detailed but often inhomogeneous description of some quantities to a rougher 
but more homogeneous representation, by averaging out inhomogeneities via 
homogenization techniques (see e.g., [5]). 

It is worth pointing out that the dichotomy fine vs. coarse scale does not 
necessarily imply a parallel dichotomy ODE vs. PDE modeling. In other words, 
it is possible that the underlying mathematical models pertain to the contin- 
uum theory at both scales (like in most of the examples recalled above) or that 
the multiscale coupling between a discrete and a continuous model is realized 
only at an approximate computational level by averaging and sampling. Con- 
versely, in the multiscale approach we propose here, the microscopic scale is 
actually a discrete one which complements the continuous flow with granular- 
ity. The resulting model is then a coupled microscopic-macroscopic one, and 
computational schemes are derived accordingly. 

The paper is organized as follows. Section[2]introduces the measure-theoretic 
framework and models pedestrian kinematics. Section [3] details the multiscale 
approach, addressing in particular the choice of microscopic and macroscopic 
parameters and their scaling. Section [4] introduces and qualitatively analyzes 
a discrete-in-time counterpart of the multiscale model. Section [5] proposes a 
numerical approximation of the equations, with special emphasis on the dis- 
cretization in space of the macroscopic scale, and explains in detail the resulting 
numerical algorithm. Section [6] discusses the results of numerical simulations 
in some case studies aimed at checking the effects of the multiscale coupling on 
the crowd dynamics predicted by the model. Section [7] finally draws conclusions 
and briefly sketches research perspectives. 

2 Mathematical modeling by time-evolving mea- 
sures 

From the mathematical point of view the mass of a d-dimensional system (d = 
1, 2, 3 for physical purposes) at time t is a Radon positive measure lit, that we 
assume to be defined on the Borel cr-algebra B(Mr). For any E G B(M. d ) the 
number fit{E) > gives the mass of pedestrians contained in E at time t > 0. 
In principle, the only further property satisfied by fi t is the c-additivity, directly 
translating the principle of additivity of the mass. 

Let T > denote a certain final time. Following [J], the conservation of the 
mass transported by a velocity field v = v(t, x) : [0, T] x R d — » R d is expressed 
by the equation 



+ V • fay) = 0, (x,t)eR d x(Q,T}, (1) 
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along with some given initial distribution of mass no (initial condition) . Deriva- 
tives appearing in Eq. ([I]) are meant in the functional sense of measures. Specif- 
ically, for every smooth test function <f> with compact support, i.e., (p G Cd(R rf ), 
and for a.e. t 6 [0, T], it results 



<j>(x) dfi t (x) = / v(t, x) ■ V(j){x) dfj, t (x), (2) 



d 
dt 



where integration-by-parts has been used at the right-hand side. A sufficient 
condition for Eq. ([2| to be well-defined is that v(t, •) is integrable w.r.t. p t for 
a.e. t e [0, T]. 

A family of time-evolving measures {/i*}t>o is said to be a (weak) solution 
to Eq. ([I]) if, for all (j) € Co(R d ), the mapping t h-> J Rd (j){x) dp t (x) is absolutely 
continuous and satisfies Eq. |2]). In particular, the latter statement means 



<p(x)d(Xt 2 (x) - J (t>(x)dfj, tl (x) = J J v(t, x) ■V4>(x)dp t (x)dt (3) 

R d t x Rd 

for all t 1 ,t 2 £ [0, T], t 1 < t 2 , and all G C^(R d ). 



Modeling the interactions among pedestrians 

Equation ([T]) provides the evolution of the measure /it as long as the velocity is 
specified. In our case, given the absence of a balance of linear momentum, this 
implies modeling directly the field v. For this reason, our approach will result 
in a first-order model. 

First-order models are quite common in the literature, especially at the 
macroscopic scale. The velocity can be either specified as a known function 
[2"o] or linked to the density of pedestrians by means of empirical fundamental 
relations v = v(p) [l>J [T^l E2] ■ Sometimes a functional dependence on the den- 
sity gradient is envisaged, in order to model the sensitivity of pedestrians to 
the variations of the surrounding density field U\ ■ Microscopic models focus 
instead more closely on the interactions among pedestrians, normally expressing 
them in terms of generalized forces. They resort therefore to a classical New- 
tonian paradigm, in which the acceleration is modeled explicitly [T^l Q3] . We 
remark, however, that in [27l [28] the authors adopt a kinematic modeling of the 
interactions in the frame of a microscopic model. 

With the aim of setting up a model based on the mass conservation only, but 
in which the microscopic granularity complements the macroscopic dynamics, 
we cannot entirely resort either to generalized forces or to fundamental relations. 
Taking advantage of the mass conservation equation in the form ^ , which does 
not assume a priori any modeling scale, our approach will be at the same time 
kinematic, macroscopic, and focused on the strategy developed by pedestrians 
at the microscopic scale. 

To be more specific, let the velocity be expressed in the following form: 

v(t, x) := v[n t ](x) = v dcs (x) + v[p, t ](x), (4) 
the square brackets denoting functional dependence on the measure pt- 
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The function v dcs : R d -> E d is the desired velocity, i.e., the velocity that 
pedestrians would set to reach their destination if they did not experience mu- 
tual interactions. In the simplest case it is a constant field, whereas in more 
complicated situations it accounts for the presence of possible obstacles to be 
bypassed (e.g., pedestrians walking in built environments). In our approach the 
desired velocity is deduced a priori from the geometry of the domain, meaning 
that it is totally independent of the measure fj,f In other words, it can be re- 
garded to all purposes as a datum of the problem. It is not restrictive to assume 
that it has constant modulus: 

I»d«(a0| = V, Vx e R d , (5) 

where V represents some characteristic speed of the walkers. We refer the reader 
to |30j for a possible method to construct v^es- 

The function v[Ht] : M. d -> K d is the interaction velocity, that is, the correc- 
tion that pedestrians make to their desired velocity in consequence of the in- 
teractions. The non-locality of the interactions is introduced in this framework 
by deriving u[nt\ from a synthesis of the information on the crowd distribution 
around each pedestrian. Specifically, we assume 

/y — *^ 
f(\y-x\)g(a x y)' rdnt(y), (6) 
\y-x\ 

R d \{x} 

where: 

• / : R+ — > R is a function with compact support describing how the walker 
in x interacts with her/his neighbors on the basis of their distance. If 
supp / = [0, R] for some R > 0, then a neighborhood of interaction is 
defined for the point x coinciding with the ball Br{x) C E. d centered in x 
with radius R; 

• a xy G [— 7r, 7r] is the angle between the vectors y — x and Vd es (x), that is, 
the angle under which a point y is seen from x with respect to the desired 
direction of motion; 

• g : [— 7r, 7r] — > [0, 1] is a function which reproduces the angular focus of 
the walker in x. 

Integration w.r.t. \i t accounts for the mass that the walkers see, considering 
that two fundamental attitudes characterize pedestrian behavior: 

• repulsion, i.e., the tendency to avoid collisions and crowded areas; 

• attraction, i.e., the tendency, under some circumstances, to not lose the 
contact with other group mates (e.g., groups of tourists in guided tours, 
groups of people sharing specific relationships such as families or parties) . 

Focusing on one of the simplest choices, nonetheless physiologically sound, 
we suggest for / the following expression: 



f( s ) = -X[o,R r ](s) + F a sx [0 ,R a ]{s), 



(7) 
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where F r , F a > are repulsion and attraction strengths, and R r , R a > are 
repulsion and attraction radii. This form of / translates the basic idea that 
repulsion and attraction are inversely and directly proportional, respectively, to 
the distance separating the interacting pedestrians. 

As pointed out in the Introduction, interactions can be either metric or 
topological. An interaction is metric if the corresponding radius is fixed, so 
that each walker interacts with all other pedestrians within that given maximum 
distance. Conversely, an interaction is topological if the corresponding radius 
is adjusted dynamically by each walker, in such a way that the neighborhood 
of interaction encompasses a predefined mass of other pedestrians s/he feels 
comfortable to interact with. In this paper we will be mainly concerned with 
metric interactions, for both repulsion and attraction. The interested reader is 
referred to [TJ [TU] , and references therein, for a detailed discussion of metric and 
topological effects, also by means of examples and numerical simulations. 

The function g carries the anisotropy of the interactions, which essentially 
consists in that pedestrians cannot see all around them and they are not equally 
sensitive to external stimuli coming from different directions. If a € [0, 7r] is the 
maximum sensitivity angular width, a very simple form of g is 

g(s) = X{\s\<a}(s), S G [— 7T, %]. (8) 

By mollifying this function it is possible to account for the visual fading that 
usually occurs laterally in the visual field when approaching the maximum an- 
gular widttQ 

3 The multiscale approach 

The framework presented in Section [2] is suitable to obtain, as particular cases, 
models at both the microscopic and the macroscopic scale. In this section we 
first briefly review the methodology for their individual derivation, already pro- 
posed in |10j to study microscopic and macroscopic self-organization in animal 
groups and crowds. Then, exploiting the tools offered by the measure-theoretic 
setting, we merge these concepts into a unique multiscale model, in which the 
microscopic and the macroscopic dynamics coexist. 

3.1 Microscopic models 

Let us consider a population of N pedestrians, whose positions at time t are 
denoted {Pj(t)}f =1 . In this case the mass of a set E G B(R d ) is the number of 
pedestrians contained in E, that is: 

m{E) =caxd{Pj(jt) € E}, 

hence /Xt is the counting measure. We represent it as a sum of Dirac masses, 
each centered in one of the -Fj's: 

N 

IH = X) ■ 

1 In order to differentiate also the maximum angular widths of repulsion and attraction 
[5], one may generalize Eq. (j6| by replacing the product of / and g with a function of two 
variables accounting simultaneously for \y — x\ and a xy . 
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Plugging this in Eq. ^ gives 

, N N 

-^^(P J (t)) = ^ U (t,F J (i)).V^(i)), \/</>eC^R d ), (9) 

3=1 3=1 

whence, taking the time derivative at the left-hand side and rearranging the 
terms, 

JV 
3=1 

the dot over Pj standing for derivative w.r.t. t. The arbitrariness of <j> implies 

P 3 {t) = v[m]{p 3 {t)), j = i,...,N, (10) 

where we have set v(t, Pj[t)) — v[nt](Pj(t)) according to Eq. Q, therefore the 
microscopic model specializes in a dynamical system of N coupled ODEs for 
the Pj's. The coupling is realized by the measure pt in the velocity field. In 
particular, the microscopic counterpart of Eq. ^ reads 

v\j*](Pj)= E f^-P 3 \)9M ^'_ P p V 

k=l,...,N 1 fe j > 

where ajy £ I -71 ") is shorthand for the angle formed by the vectors Pk — Pj 
and Vdes(Pj)- We point out that, with the function / given by Eq. Q, the 
statement Pk ^ Pj in the above formula can be converted into the milder one 
k =/= j. Indeed one can prove that if the Pj's are initially all distinct they remain 
distinct at all successive times t > (see [9] for technical details). 

3.2 Macroscopic models 

Macroscopic models are based on the assumption that the matter is continuous, 
thus the measure pt is absolutely continuous w.r.t. the <i-dimensional Lebesgue 
measure C d , p t < C d . Radon-Nikodym's Theor em asserts that there exists a 
function p(t, ■) £ L 1 1 oc (E rf ) such that 

dp t = p(t, -)dC d , p(t, ■) > a.e., (11) 

called the density of pt w.r.t. C d . In our context p(t, x) represents the density 
of pedestrians at time t in the point x. 

Using p, the mass conservation equation ^ rewrites as 

~ J p{t, x)<f>(x) dx = J p(t, x)v(t, x) ■ V0(x) dx, G C£(M d ), (12) 

R d R d 

namely a weak form of the continuity equation 

g + v-(H = o. (13) 

The interaction velocity specializes as 

f(\y - x\)g(a xy )- -p(t, y) dy 

\y-A 

where it should be noticed that the domain of integration may now indifferently 
include or not the point x because {x} is a Lebesgue-negligible set. 
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3.3 Multiscale models 

If the measure p, t is neither purely atomic nor entirely absolutely continuous 
w.r.t. C but includes both parts, we get models that incorporate the micro- 
scopic granularity of pedestrians in the macroscopic description of the crowd 
flow. More specifically, we consider 

IH = Bmt + (1 - 9)M t , (14) 

where 

N 

m t = 2J S Pj(t), dM t (x) = p(t, x) dx 
3=1 

are the microscopic and the macroscopic mass, respectively. The parameter 
9 G [0, 1] weights the coupling between the two scales, from 9 = corresponding 
to a purely macroscopic model to 9 = 1 corresponding to a purely microscopic 



model. In Eq. (14) no scaling parameters explicitly appear, but we anticipate 



that they will arise naturally from our next dimensional analysis (cf. Section 



3.4) 



Using the measure (14), the mass conservation equation |2]) takes the form 
of a mix of microscopic and macroscopic contributions: 

p(t, x)cf){x) dx = 



d 
dt 



(05>(p,-(*)) + (i-f) / / 

N 

622v(t, Pj(t))- V^{Pj(t)) + (l-e) / p(t, x)v(t, x) ■ V<t>(x) dx, V^eCo 1 ^), 



j'=i 



formally a convex linear combination of Eqs. ([9]), ( 12 ). The interaction velocity 
v [fit] is now given by 

v\ 1H \{x) = e J2 f(\Pk(t)-x\)g(a xPk{t) ) Pk ^l X 

k=l,...,N ^ k[l) ' T| 



1 y — *^ 

+ ( 1 -9) / f(\y-x\)g(a xy )- p{t, y)dy, 

\y - x\ 



therefore it coincides neither with the fully microscopic nor with the fully macro- 
scopic one. This definitely makes the overall dynamics not a simple superposi- 
tion of the individual microscopic and macroscopic dynamics. 

It is worth noticing that the point x may or may not be one of the po- 
sitions of the microscopic pedestrians. Computing v[pt] for x = Pj(t) shows 
that the interaction velocity of the j-th pedestrian accounts not only for other 
microscopic pedestrians contained in the neighborhood of interaction but also 
for the macroscopic density distributed therein, which represents some crowd 
whose subjects are not individually modeled. Specifically, the term responsible 
for this is 

f(\y- Pj {t)\)g{a P . {t)y ) ^^- [ P {t, y)dy, (15) 
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that we may regard as the macroscopic contribution to the microscopic dynam- 
ics. Analogously, computing v[pt\ for x different from all of the Pj's shows 
that the interaction velocity of an infinitesimal reference volume centered in x 
depends not only on the density distributed in the neighborhood of interaction 
but also on the microscopic pedestrians therein, which play the role of singular- 
ities in the average crowd distribution due to the granularity of the flow. The 
corresponding term is 

E f(\ P k(t)-x\)g(a xPk{t) ) ^- X (16) 

k=l,...,N l ^ k[T> Xl 

which gives the microscopic contribution to the macroscopic dynamics. 



3.4 Dimensional analysis 

In order to scale correctly the microscopic and the macroscopic contributions, 
it is convenient to refer to the non-dimensional form of the model. For this, let 
us preliminarily notice that the main quantities involved in the equations have 
the following dimensions: 

• [t] = time 

• [x] = length 

• [vdes] = [v] = length/time 

• [/] = lcngth/(time x pedestrians) 

• [/it] = pedestrians 

• [p] = pedestrians/length d 

where "pedestrians" is actually a dimensionless unit. Additionally, g and 9 are 
dimensionless. Let L, V, g be characteristic values of length, speed, and density 
(in particular, V may be the desired speed introduced in Eq. ^) to be used to 
define the following non-dimensional variables and functions: 

x V 1 

f*(s*) = ~f(Ls*), p*(t*,x*) = ~p(yt*,Lx*y P*(t*) = ~P j (~t*y 

Notice that, due to the choice of V as characteristic speed, the dimensionless 
desired velocity v^ es turns out to be a unit vector. 

In more detail, the non-dimensional mass measure p*, is given by 

dp* t ,(x*) = dp,L t ,(Lx*) 

= eZjhpf^)(Lx*) + (l-e)gL d p*(t*, x*)dx* 

= ej2 j s PjHt , ) (x*) + (i-e)Ap*(t*, x*)dx* 

= 6m* t , {x*) + (1 - 6) AM*, (x* ) 
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where we have set A := gL d and we have recognized the dimensionless micro- 
scopic and macroscopic masses: 

N 

m t' = % (** ) ' dM t* ( X * ) = P* (** ' ' ) dx * ■ 

We notice that the coefficient A has unit [A] = pedestrians, therefore it is a non- 
dimensional number fixing the scaling between the microscopic and the macro- 
scopic mass. It says how many pedestrians are represented, in average, by a 
unit density p* in the infinitesimal reference volume dx* . 

Remark. The measure 



M*« = 0rn* t , + (1 - 0)AM t * (17) 

can be read as a linear interpolation between the microscopic and the macro- 
scopic mass via the parameter 8, provided m*, (M. d ), M t *„ are, up to scaling, 
the same mass, i.e., m*, (M. d ) = AM t *. (R d ). As we will see later (cf. Corollary [|]), 
in the multiscale model the microscopic and macroscopic mass are individually 
conserved in time, hence this can be achieved by setting 

_ m* (R d ) N 

M*(R d ) M*{R d ) 1 ; 

as long as < N, M^(R d ) < +oo. 

In the following we will invariably refer to the non-dimensional form of the 
equations, omitting the asterisks on the non-dimensional variables for brevity. 



4 Discrete-in-time model 

In this section we derive a discrete-in-time counterpart of the multiscale model, 
that will help us gain some insights into the qualitative properties of the math- 
ematical structures previously outlined. In addition, it will serve as a first step 
to devise a numerical scheme for the approximate solution of the equations. 

Let At n > be a possibly adaptive time step and let us introduce a sequence 
of discrete times {i n } n >o such that i = and t n+ i — t n = At n . Denoting 
(i n := fit„, from Eq. ^ with the choice t\ = t n , ^ = t n +x we get 

t„+i 

4>{x) dp n+1 (x) - J 4>(x) dp, n {x) — J j v(t, x) ■ V 4>{x) dp t {x) dt 



= At n J v(t n , x) ■ V(f>(x) dfi n (x) + o(At n ), 

whence 

4>(x) djj, n+ i{x) = / [(j)(x) + At n v(t n , x) ■ V0(x)] dfi n (x) + o(At n ). 
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At this point let us explicitly assume that p n (M. d ) < +00. If v(t n , •) is p n - 
uniformly bounded then 4>(x) + At n v(t n , x) ■ V(f>(x) = 4>{x + At n v(t n , x)) + 
o(At n ), thus 

4>{x) dfJLn + i(x) = / 4>(x + At n v{t n , x))d^ n {x) +o(At n ). 



Defining the flow map ^y n {x) := x + v(t ni x)At n and neglecting the term o(At„), 
we are finally left with 



(x)dfi n+ i(x) = j c/)(j n (x))dfj, n (x), (19) 

R d R d 

which makes sense actually for every bounded and Borel function 4>. Choosing 
4> = Xe for some measurable set E e B(R d ) entails 

M„+i(£) = Vnh-\E)), VE e 6(M d ), 

meaning that /x n +i is the push forward of /Lt n via the flow map 7„, also written 



Mn+i = 7n#Mn- Equation (19) provides a discrete-in-time counterpart of Eq. 
([3]) . Obviously, it requires to be supplemented by an initial condition /iq in order 
for the sequence {n n }n>i to be recursively generated. 

Notice that, with the velocity field Q, it results v(t n , x) — v\p, n ]{x) with in 
particular: 



k=l,...,N 



+ (1-9)A f f(\y-x\)g(a xy )f^- Pn (y)dy, (20) 

J \y-x\ 

where PJ} := P k (t n ) and p n := p(t n , ■). 

Preserving the multiscale structure of the measure 

Recall that in the multiscale model we assumed that our measure is composed 
by a microscopic granular and a macroscopic continuous mas^J Of course, this 
is just a formal assumption made to write the model. From the analytical point 
of view, it need be proved that such a measure can be actually a solution to our 
equations. 



Set m n :— m tn , M n := M tn , so that, owing to Eq. (17 1, the measure p n can 
be given the form 

fx n = 6m n + (1 — 6)AM n . (21) 

The following result clarifies the role played by the flow map j n in preserving 
the multiscale structure of p n after one time step. 



2 With respect to the general structure of a measure as provided by Ricsz's Theorem, this 
means that we are in particular excluding the Cantor's part. 
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Theorem 1. Let a constant C n > exist such that 

C d (i-\E)) < C n C d (E), V£ g B(R d ). (22) 



Given p n as in Eq. (21), there exist both a unique atomic measure m„ + i and 
a unique Lebesgue- absolutely continuous measure M n+ i, with a.e. nonnegative 
density, such that /J, n +i = Om n+ i + (1 — 9)KM n+ \. 

Proof. 1. Using the linearity of the operator 7 n #-, the measure fi n +i is given 

by 

Mn+i = 0{-Yn#rn n ) + (1 - 6»)A( 7 „#M„). 

2. Let us define 

m n+ i := 7„#m„. 

A direct calculation shows that such a m n +\ is in fact an atomic measure. For 
any measurable set E G S(R. ) we compute: 



N 

[ ln #m n ){E) = m^-^E)) = ^ S P? (E)) 

3=1 



= card{ 7 „(P J ") G £} 

JV 
3=1 

hence m n+ i as defined above is in turn a combination of Dirac masses centered 
in the new positions {P™ +1 }jLi given by 

P; +1 := 7 „(P; 1 ) = Pp + v{p n ](P?)At n . (23) 

3. Analogously, let us define 

M n+1 := 7n #M„. 

We claim that, under the hypothesis of the theorem, this measure is absolutely 
continuous w.r.t. C d . To see this, let E G B(R d ) be such that C d {E) = 0. Then 
^ d (7r7 1 (^)) < C n C d (E) — 0, whence, using that M n <C £ d by assumption, we 
get M n+1 (E) = M n (^ 1 (E)) = and the claim follows. 

4. To show the non-negativity of the density of M n+ \ we take the Radon- 
Nikodym derivative. Then we discover, for a.e. x, 

M n+1 (B r (x)) 1 f 
p n+ i{x)= hm = hm ^ / p„(y)dy, 

■r^0+ £ a (±S r {X)) r->0+ W<jr a J 

where is the volume of the unit ball in K d . Since p„ is a.e. nonnegative by 
assumption, the same holds for p n +\ and we are done. 

5. Finally, uniqueness of m n+ \, M n +\ is implied by the uniqueness of the 
Radon-Nikodym decomposition of a measure. □ 

The proof of Theorem [T] is constructive, indeed it shows explicitly how to 
obtain the measure p n +i starting from p n : one simply pushes forward separately 
the microscopic and the macroscopic mass via the common flow map 7 „. 

By referring to the results proved in 29 , we can state some additional 
properties of the measure p n . 



MULTISCALE MODELING OF GRANULAR FLOWS 



13 



Corollary 2. If 7„ satisfies (22) /or all n > and the initial measure /xq 



complies with the form (21) i/ien 



i. there exists a unique sequence of atomic measures {m n } n >i and a unique 
sequence of positive Lebesgue- absolutely continuous measures {M„}„>i 



such that fx n has the form (21 ) for all n > 1; 

2. the measure fi n satisfies the following conservation law: 

» n+ i(E) - fi n (E) = - [fi n (E \ 7~ 1 (i?)) — ti n (^\E) \ E)] (24) 

for all E € B(M. d ). In particular, both m n and M n satisfy this law sepa- 
rately at each time step; 

3. if po € L^ c (R d ) then p n € L^ c (R d ) for all n > 1, with moreover 

n-i 

esssup |p„(x)| < TT Cj ess sup \po(x)\ 

x£E j =0 x£E 

for all compact set E C K. , where the Cj 's are those appearing in the 
statement of Theorem^ 

Proof. 1. Existence and uniqueness of microscopic and macroscopic masses have 
been proved in Theorem[l]for one time step, hence they follow for all times n > 1 
by induction. 

2. In view of the er-additivity of the measure we have, for all E € B(R d ), 

Mn+l(£) = Vnhn'iE) H E) + ^~\E) \ E). 

Subtracting p, n (E) at both sides and collecting conveniently gives 

Mn+l(£) ~ »n(E) = - [l*n(E) - ^{^{E) fl E)] + ^-'(E) \ E) 

whence, observing that E — (^^(E) n E) U (E \ 7 r ^ 1 (-E)) with disjoint union, 
the thesis follows. Since we know from Theorem [T] that m„ +1 , M„ +1 are in turn 
generated by push forward with j n , this reasoning can be repeated to find that 
each of them fulfills the very same conservation law. 

3. In [29] it is proved that p n e L^ c {R d ) implies p n+1 € L^ c {R d ) as well, 
with esssup^g^ |p„+i(a;)| < C n esssup^g^ |/? n (a;)|. Thus proceeding by induc- 
tion from n = we get the result. □ 

Some comments on the results of this section are in order. 

(i) The main assumption of both Theorem [l] and Corollary [2] is that the flow 



map 7„ satisfies Eq. (22 1. In general it may be hard to check the validity 
of this property directly but in |29j it is proved that a sufficient condition 
for it to hold true is that the velocity v[p, n ] be Lipschitz continuous and 
that the time step be chosen so that At n Lip(v[/i n ]) < 1, n = 0, 1, 2, 



(ii) Equation ( 24 ) states that the variation of the mass of a set E in one time 
step is due to the net mass inflow or outflow across dE. Indeed, E\j^ 1 (E) 
is the subset of E which is not mapped into E by 7„ (outgoing flux) , and 
7~ 1 (i?) \ E is the subset of R d \E which is mapped into E by j n (ingoing 
flux). 
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(iii) Choosing E = R d in Eq. ([24} yields p n+ i(R d ) = fi n (R d ), i.e., the conser- 
vation of the total mass at each time step. If /io(M. d ) < +00 then the mass 
is finite for all n, therefore, up to normalization, the /Lt„'s may be regarded 
as probability measures. Furthermore, since Eq. (24 1 applies separately 
also to m n and M n , both the total microscopic and the total macroscopic 
masses are conserved in time. 



5 Numerical approximation of the equations 

As anticipated at the beginning of Section[4j the discrete-in-time model provides 
a first discretization of the equations, which would be sufficient for tracking the 



microscopic mass (cf. Eq. (231). However, the macroscopic mass requires a 



further discretization in space in order to come to a full approximation of the 



density p. Notice that, in principle, one may refer to Eq. (13) and rely on the 
wide literature on numerical methods for nonlinear hyperbolic conservation laws 
[24j . Nevertheless, aside from the intrinsic complication due to the multidimen- 
sional nature of the equations, this strategy poses several nontrivial technical 
difficulties. For example, it demands a correct definition of the convection ve- 
locity (i.e., formally the derivative of the flux pv with respect to the density) 
in presence of nonlocal multiscale fluxes, as well as a consistent formulation of 
entropy-like criteria for picking up physically significant solutions. All these 
issues are instead bypassed if one maintains the measure-theoretic formalism. 

For the discretization in space of the density p n we partition the domain 
in pairwise disjoint d-dimensional cells Ei E B(M. d ), where i £ 1 d is an integer 
multi- index, sharing a characteristic size h > such that C d (Ei) — > for all i 
when h — > + (for instance, h ~ diamE^). Every cell is further identified by 
one of its points Xi, e.g., its center in case of regular cells. 

We approximate p n by a piecewise constant function p n on the numerical 
grid: 

p n (x) = p™, Vi £ £i, 

where p™ > is the value that p n takes in the cell Ei. Consequently, the measure 
M n is approximated by the piecewise constant measure dM n — p n d£ d , which 
entails the approximation p n — 0m n + (1 — 6)AM n for p, n . 

Analogously, we approximate the velocity v[p n ] by a piecewise constant field 

v[p n ](x) = v™, Vie E u 

where the values v" G K d are computed as vf = v[p n ](xi). The discretization 
of the velocity gives rise to the following discrete flow map: 

%(x) = x + v[p n ](x)At n , 

which turns out to be a piecewise translation because t>[/2„] is constant in each 
cell. 

Finally, we look for a piecewise constant approximation M n+ \ of M n+ \ by 
imposing the push forward of M n via the flow map 7„: 



M n+1 (E) = M n (% 1 (E)) 1 V£ e B(R d ). 
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In particular, choosing E — Ei yields 

P7+1 = £^E~) £^ d (^n7n(£ fc )), V^Z", (25) 

which provides a time-explicit scheme to compute the coefficients of the density 
p n +i from those of p n . Notice in particular that j n (Ek) is simply the set Ek + 
v%At n . 

Notice that this scheme is positivity-preserving, in the sense that p n > 
implies p n +i > as well, hence, by induction, po > implies p n > for all 
n > 0. Such a basic property is not as straightforward in usual numerical 
schemes for hyperbolic conservation laws. Indeed, unless suitable corrections 
are implemented, the latter may develop oscillations leading to locally negative 
approximate solutions even when the exact solution is not expected to be so. 

Furthermore, considering that j n is a translation in each grid cell and using 
the invariance of the Lebesgue measure under rigid transformations, we deduce 

f p n+1 (x) dx=J2 Pk E Cd &n HEi) n E k ) = J2 Pk&iEk) = [ Pn{x) dx, 

thus the approximate macroscopic mass M n is conserved in time. 

The quality of the spatial discretization described above with respect to the 
refinement of the grid, in the case of regular flow maps, is provided by the 
following result. 

Remark. At this point we assume explicitly that the domain of the problem is a 
bounded set 51 C K d , which for all fixed h > is partitioned with a, finite number 
of grid cells (however tending to infinity when h — > + ). The multi-index i of 
the grid cells runs in a finite subset I C Z d . 

Theorem 3. Assume that 7„ is a diffeomorphism and let h, At n be sufficiently 
small and satisfying 

max^K|<l. (26) 

Then: 

(i) There exists a constant C n > 0, independent of h, such that 

\M n+1 (Ei) - M n+ i(E l )\ <C n i J \ Pn {x) - p n (x)\dx + h 

(ii) Ifv[p n ](x) is uniformly bounded, there exists a constant C' n > 0, indepen- 
dent of h, such that 

max\M n {Ei) - M n (E t )\ < max |M (£i) - M (Ei)\ + C' n h d . 

Proof. See [12]. □ 

In order to gain some control over the error introduced by the spatial dis- 
cretization, Theorem [3] requires the CFL-like condition ( 26 1 to be satisfied at 



each time step, similarly to numerical schemes for hyperbolic conservation laws. 
However there is a remarkable difference from their CFL condition, namely that 



Eq. ( 26 ) involves directly the flux velocity and not the convection velocity. 



MULTISCALE MODELING OF GRANULAR FLOWS 



16 



The algorithm 

Here we detail the numerical algorithm stemming from the above scheme, that 
we use for simulations. 

The algorithm combines a microscopic and a macroscopic part. The former 
handles the evolution of pedestrian positions, updating a vector which stores the 
values Pf € R d . The latter manages instead the evolution of the density, and at 
every time step it updates the values p™ at the grid cells. The two models evolve 
by means of the same velocity field 5[/Xn]j thus guaranteeing coherence of the 
final solution. This conceptual scheme is motivated by Theorem[T] The velocity 
field must be defined at pedestrian positions {-P"}jLi f° r the microscopic part 
and at the grid cells {^ijigi for the macroscopic part. 

Let us introduce the following superscripts: 

• micro: quantities defined at pedestrian positions; 

• macro: quantities defined at grid cells; 

• micro-for-micro: microscopic quantities computed at pedestrian positions; 

• micro-for-macro: microscopic quantities computed at grid cells; 

• macro-for-micro: macroscopic quantities computed at pedestrian posi- 
tions; 

• macro-for-macro: macroscopic quantities computed at grid cells. 
The algorithm consists of the following steps. 

1. Initialization. We fix the number N of microscopic pedestrians that we 
want to model, we define their positions, and we compute the coefficients 
pi of the initial density according to a local average of the microscopic 



mass, taking the scaling (18) into account. More precisely, we set: 



_ m (B c (x t )) . 
~~ AC d (B^(x l )) ' ieZ > 

where mo is the microscopic mass at the initial time and B^(xi) is the ball 
centered in the center of the grid cell Ei with radius £ > 0. The latter is 
tuned depending on the positions of the microscopic pedestrians, in such 
a way that the relation A = m (R d ) / AI (R d ) be satisfactorily fulfilled in 
the numerical sens^] (Mq being the approximate macroscopic mass at the 
initial time). 

2. Microscopic part. At time t = t n we compute the sum at the right-hand 
side of Eq. ( 20 ) for x — obtaining 



-micro-for-micro . = ^^(p") 



3 

The same computation performed for x — Xi gives instead 

^micro-for-macro := u[m n ]{Xi), 



cf. Eq. (16), which will be shared with the macroscopic part of the code. 



3 Notice that if one replaces B^{x{) with the cell Ei then the measures mo, Mo satisfy the 
scaling (JTsJ) exactly. However, averaging on a neighborhood a bit larger than a single grid 
cell is essential in order to have a macroscopic density really distributed in space rather than 
clustered in grid cells. 



MULTISCALE MODELING OF GRANULAR FLOWS 



17 



Macroscopic part. At the same time instant t = t n we numerically evaluate 



the integral at the right-hand side of Eq. (20 1 for x = Xi, using the 



approximate density p„ in place of p n . This way we obtain 

-macro-for-macro ;= ^M n ]{ Xi ). 

Next we compute the same integral for x = P" , which yields 

-macro-for-micro . = ^M n ](P^), 



cf. Eq. (15). In particular, the integrals involved in £ macro - for - macr ° an( j 
~macro-for-micro are numerically evaluated via a first order quadrature for- 
mula. This component of the velocity field will be shared with the micro- 
scopic part of the code. 

4. Desired velocity. If the velocity field Wd cs is given analytically, the com- 
putation of w™ s cro := Vdcs(Pj l ) and of f^" := Wdcs(^i) is immediate. If 
instead «dcs is defined on the numerical grid only, for instance because it 
comes from the numerical solution of other equations [30], then w™" is 
computed by interpolation. Since we are assuming that all macroscopic 
quantities are piecewise constant, we coherently choose a zeroth order 
interpolation. 

5. Overall velocity. We assemble the previous pieces as 

« micr ° : = i\fin](If) 

micro , n ~ micro- for- micro 1/1 n\ \ -macro-for-micro 

and analogously 

S macro : = 5[£„](:e<) 

n , macro , r\ —micro- for-macro 1/1 n\ a -macro- for- macro 
= v dcs + I 1 - V)l\v 

6. Computation of At. We compute the largest time step At allowed by 



condition ( 26 ) for the macroscopic velocity field v 1 ' 



7. Advancing in time. We update pedestrian positions and density according 



to Eqs. (23), (25) by means of w mlcro and f) macl °, respectively. 



Remark. No matter what the value of is, the two approaches always coexist. 
If 9 = the macroscopic scale is leading, and the microscopic pedestrians are 
simply driven by the macroscopic velocity field. This is the classical way to 
see flowing (Lagrangian) particles in a fluid, whose motion was previously com- 
puted. Conversely, if 6 = 1 the microscopic scale is leading, and the evolution of 
the macroscopic density is reliable only if the number of microscopic pedestrians 
is sufficiently large. 



6 Numerical tests 



In this section we present the results of numerical simulations performed with the 
model and the algorithm described above. As natural for pedestrian flows, we 
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Table 1: Summary of the parameters used in the numerical tests 



Test N. 


9 


N 


A 


F r 


F a 


R r 


Ra 


1 


[0, 1] 


100 


10 


0.1 





0.5 


N/A 


2 


[0, 1] 


10, 100 


10, 100 


0.1 





0.25 


N/A 


3 


0, 0.3, 1 


30 


30 


0.1 





see text 


N/A 


4 


0.3 


25 + 1 


80 


0.05 


0.4 


1.5 


1.5 



deal with two-dimensional (d = 2) bounded domains, say f2 C R 2 , confining the 
attention to the restriction measure /iti_f2. This means that the mass possibly 
flowing out of the domain is considered as lost, i.e., it no longer affects the 
computation. 

Sometimes we will deal with domains with obstacles understood as internal 
holes. They require a careful handling of the velocity at their boundaries so 
as to prevent it from pointing inward (which would imply unrealistic outflow 
of mass). In order to have the mass bypass the obstacles, the velocity Q is 
projected onto a space of admissible velocities V a dm, which can be defined in 
several ways depending on the pedestrian behavior one wants to model. Our 
choice for the next examples is 

Vadm — {v €M. d : v ■ n > at every obstacle boundary}, 

where n is the outward normal unit vector at the obstacle boundaries. This 
corresponds to setting to zero the normal component of the velocity Q in case 
it points into an obstacle. A different possibility is to set to zero both the normal 
and the tangential component if the first one points into an obstacle. In the 
former case pedestrians can slide along the obstacle walls following the tangential 
velocity, whereas in the latter case they remain still against the obstacles until 
no longer pushed by flowing neighbors. This choice may model, for instance, 
a more relaxed condition in which walkers are not in a hurry to reach their 
destination. 

Concerning the parameters, we assume a = n/2 (frontal interaction) in all 
tests, which is suitable for the most common situations encountered in pedes- 
trian flow. We also assume no attraction between group mates but in the last 
test (Test 4) in which we model the dynamics of a group of people following 
a leader. Table [T] summarizes the values of all other parameters used in the 
numerical tests. 



Test 1: Dynamics of the interactions 

In this first test we study the effect of the multiscale coupling on the rearrange- 
ment of a crowd subject only to internal repulsion. The goal is to show that 
it is possible to obtain a perfect correspondence between the microscopic and 
the macroscopic dynamics in some simple cases, which originally motivated and 
justified the possibility of a coupled multiscale approach [TU] . 

To this purpose we switch the desired velocity off, so that the velocity v[/j,t] 
coincides with the interaction velocity z-4/it]. (Actually, in order to compute the 
angle a xy in Eq. ( 20 ) we conventionally assume Vdcs to be constantly directed 
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1 2 3 4 5 6 

(c) 

Figure 1: Test 1. (a) Initial condition, 
(b) the purely macroscopic model, (c) 
the multiscale model with 8 = 0.3 
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(b) 
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1 2 3 4 5 6 

(d) 

Crowd distribution at time T = 1 with 
the purely microscopic model, and (d) 



along the horizontal axis, so as to define what is ahead and what is behind). 
Pedestrians are initially arranged in a square-shaped equally-spaced formation, 
see Fig. [IJl Due to the frontal repulsion, we expect the frontal part of the 
group to stand still and the rear part to stand back from the group mates 
ahead. We compare the expansion dynamics of the group as predicted by the 
macroscopic (9 — 0), the microscopic (8 — 1), and the multiscale (9 = 0.3) 
model. The simulation runs until the final time T = 1 is reached. Notice that 
the configuration assumed at that time is not an equilibrium of the system. 
Results are shown in Figs. [lp-d. 

The main features of the dynamics outlined above are caught at all scales. 
In particular, the effect of the only frontal repulsion is visible at the head of 
the group, where pedestrians stay aligned on a vertical line as they are initially 
because there is none in front of them. This clearly shows up looking both at 
the density distribution at the macroscopic scale (Fig. [T]d) and at the individual 
pedestrians at the microscopic scale (Fig. [lj). 

Of course, this does not mean that either scale has no influence at all on 
the other. For instance, as an interesting effect of the microscopic scale driving 
the macroscopic dynamics, we notice some kind of "density holes" near every 
microscopic pedestrian in the limit of the purely microscopic model (Fig. [lp). 
They are actually small areas of very low density, caused by that microscopic 
repulsion has a great impact at the macroscopic scale. Recall indeed that the 
microscopic granularity is seen as a singularity in the average crowd distribution, 
and that for 8 = 1 the evolution of the macroscopic density is fully ruled by 
the microscopic scale. With the multiscale model (Fig. [TJi with 9 = 0.3) the 
hole effect is instead limited, and a good compromise between the two scales 
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is reached. Furthermore, in Fig. [TJi (multiscale) pedestrians are less scattered 
than in Fig. JTJi (microscopic), meaning that the contribution of the macroscopic 
scale on the overall dynamics has, in a sense, a homogenizing effect. Conversely, 
in Fig. [TJi the macroscopic density is more broken than in Fig. [IJd (macroscopic) , 
thus the microscopic scale destroys the macroscopic smoothness and introduces 
a non-negligible granular effect in the overall dynamics. 

To investigate more in depth the intercorrelation between the scales we con- 
sider now how the moments of inertia of the mass distribution depend on the 
coupling parameter 6, fixing all other parameters as indicated in Tablc[T] Indeed 
from classical mechanics it is known that moments of inertia provide quantita- 
tive information on the shape of the group. 

Let xg be the center of mass of the crowd at the final time T: 



XG = — TFa / xd VT{x), 



<> 



then we consider the following three moments of inertia around xq'- 

/f = J \(x-x G )-ifdrr{x), /£ = J |(x-z G )-j| 2 dMz), = # + 

i, j being the unit vectors in the direction of the horizontal and vertical axis, 
respectively. If and 1% refer to stretching or shrinking of the group in the 
horizontal and vertical direction, respectively, whereas I G accounts for the global 
distribution of the crowd around its center of mass. By replacing [i? in the above 
formulas with the measure mx (Mt, resp.) it is possible to study the analogous 
moments of inertia of the sole microscopic (macroscopic, resp.) mass, that we 
denote by I^ 2 G (l£* 2 Gj _resp.) 



The graphs of Fig. l2l show the trend of the functions 0^1 
I?'*"'", and 6 i-» j™>-«^ Notice t w due to Eq. rfPfb, the 



1 , - 



l g . nuwuc unau, uuc TO Eq. (17), the moments of 
inertia of the multiscale mass are linear interpolations of the corresponding 
moments of inertia of the microscopic and the macroscopic masses. The latter 
are therefore also plotted in the graphs for reference. The most relevant fact 
is that the multiscale moments of inertia are almost constant with respect to 
(aside from small border effects, especially about 8 = 1), which indicates that 




the rearrangement of the mass is basically the same at all scales. Therefore the 
microscopic and the macroscopic dynamics arising from pedestrian interactions 
are compatible with each other and make it possible a coupled approach by scale 
interpolation. 

The graphs of Fig. [i] show the trend of the mappings t H> Iq ' M ' p for the 
three values of 9 used in Fig. [I] namely 9 = (Fig. 9 = 0.3 (Fig Efr ), 
and 9 = 1 (Fig. [3J;). As pointed out in the Remark at the end of Section [5 the 
microscopic and the macroscopic scale always coexist and exchange information. 
In particular, by comparing Figs. [3^, and p3p it can be noticed how the scale 
coupling realized in the model produces coherent results at both scales even when 
the dynamics is fully ruled by either of them only. No significant qualitative and 
quantitative differences are observed in both the fully macroscopic {9 = 0, Fig. 
3ja) and the fully microscopic (9 = 1, Fig. [3p) case, meaning that there is no 
detachment of the two evolutions even when only one of them is actually the 
leading one. If this might be quite classical for the (Lagrangian) evolution of 
microscopic particles driven by a macroscopic flow, we stress that it is definitely 
by far less classical and obvious for the (Eulerian) evolution of a macroscopic 
flow driven by microscopic particles. 



Test 2: Average outflow time 

In this test we address the case of a crowd leaving a room through a door in 
normal (i.e., no panic) conditions, and we investigate the influence of the coupled 
microscopic and macroscopic effects on the estimated average outflow time. 
This will provide meaningful insights into the way in which the microscopic 
granularity works within the macroscopic flow. The scenario of the simulation 
is depicted in Fig. [4^, for the parameters listed in Table [l] The (dimensionless) 
door width is 0.5. 

Let Ct := [0, 3] x [0, 4] be the room that pedestrians are leaving. We consider 
the following average outflow time: 

T 

T» vc := — / tF(*) #> 
where J°(t) is the integral flux through the door (taken positive when outgoing) 
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3 8 6 9 

(a) (b) (c) 

Figure 4: Test 2. (a) A crowd leaving a room through a door, initial condition 
(top) and underway outflow (bottom), (b) Average outflow time as a function 
of 9 for a crowd of 10 pedestrians and (c) 100 pedestrians 



at time t. The final time T > is chosen so large that the room is defi- 
nitely empty, i.e., ht(&) = 0. Considering that J*(t) = — ^fi t (&) and using 
integration-by-parts, T£ vc can be given the following form: 

T 

T^^^Lmdt, (27) 

Mo(") J 
o 



hence it is actually an outflow time weighted by the percent mass of crowd that, 
at each time instant, still has to leave the room. 

The graphs of Figs. [4Jd, c show the trend of the function 9 h-> T£ vc for a 
small crowd of 10 pedestrians and a large crowd of 100 pedestrians. In both 
cases, the two further curves 9 i-» T™ c , 9 H> T^f e , computed by replacing fx t in 



Eq. (27 1 with either m t or M t , are plotted for suitable reference. Again, due to 



Eq. (17), the multiscale average outflow time (27 1 is the linear interpolation of 
the corresponding microscopic and macroscopic times via the parameter 9. 

The trend of the T ave 's is qualitatively similar for both the small and the large 
crowd, in particular it is decreasing with 9. This elucidates the role played by a 
more and more influential microscopic granularity within the macroscopic flow: 
the more the multiscale coupling is biased toward the microscopic scale, the more 
fluent the crowd stream becomes (and consequently the average outflow time 
decreases). This is justifiable considering that 9 can be viewed as the percent 
mass shifted from the macroscopic to the microscopic scale in consequence of 
the multiscale coupling. Subtracting interacting macroscopic mass from the 
system progressively reduces the action of the macroscopic interactions while 
enhancing that of the microscopic ones. Since the latter are less distributed 
in space, because the microscopic mass is clustered in point singularities, this 
ultimately results in fewer deviations from the desired velocity and the desired 
paths. 

Test 3: Pedestrian flow through a bottleneck 

In this test we investigate the ability of the multiscale model to reproduce several 
flow conditions occurring when two groups of pedestrians, walking toward one 
another, share a narrow passage (e.g., a door). 
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Figure 5: Test 3. Clogging at the bottleneck (left) and corresponding macro- 
scopic fluxes (right) arising with the fully macroscopic dynamics 



From the modeling point of view it is necessary to handle two interacting 
populations of walkers, which will be done via two mass measures \jP t , p = 1, 2, 
each obeying Eq. ([I]). Either population has its own desired velocity v p cs and 
interacts with the opposite population through the interaction velocity, which 
now depends on both the juf's. Specifically, denoting by p* the conjugate inde^ 
of population p, we set 

iftf, nf] = (i - e)i/>?] + e^V ], (28) 

where: 

• v v \[iP t \ is the endogenous interaction, i.e., the interaction with pedestrians 
of one's own population; 

• v pp [fi P ] is the exogenous interaction, i.e., the interaction with pedestrians 
of the opposite population; 

• € [0, 1] is a dimensionless number fixing the strength of the exogenous 
against the endogenous interaction. 

Both the endogenous and the exogenous interaction velocities are formally com- 
puted as in Eq. ([6]), except that the exogenous one is integrated with respect to 
fi P . In addition, the exogenous interaction radius R pp need not be the same 
as the endogenous one RP if interactions with opposite walkers require more 
promptness than interactions with group mate^] 

For this test we let O = 0.65, thus 65% repulsion is exogenous and 35% is 
endogenous. Repulsion radii are R p = 0.2, R pp = 0.35, to be compared with 
the unit width of the narrow passage. Other relevant parameters are listed in 
Table [TJ The setting of the problem is displayed in the snapshots of Figs. |5j|HJ 
[6j in particular the blue crowd with red microscopic pedestrians, say population 
1, walks rightward whereas the yellow one with green microscopic pedestrians, 
say population 2, walks leftward. 

Let us begin from the case 9 — 0, with the macroscopic scale leading the 
dynamics. The bottleneck tends to clog (Fig. [5^,): no density nor microscopic 

4 That is, p* = 2 if p = 1 and p* = 1 if p = 2. 

5 This simply corresponds to the function / having different supports in the expressions of 
yP and v™* . 
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Figure 6: Test 3. Alternate flows at the bottleneck in the multiscale model 
(9 = 0.3). Negative values of the density of the population walking rightward 
are for pictorial purposes only 
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Figure 7: Test 3. Multiscale model (9 = 0.3): (a) macroscopic fluxes across the 
bottleneck, (b)-(c) number of microscopic pedestrians crossing the bottleneck 
rightward and leftward, respectively 



pedestrians flow through, except for a small mass passing initially when the 
passage is still free. This is well confirmed by the time trend of the macro- 
scopic flux across the bottleneck (Fig. [5]d) : that of population 2 is permanently 
zero for t > 4.5, whereas that of population 1 oscillates between small positive 
and negative values, which implies that population 1 is pushed backward by 
population 2 as soon as it tries to cross. 

By increasing 9 to an intermediate value between and 1, a multiscale 
coupling is realized. For 9 = 0.3, the resulting dynamics is depicted in Fig. [6] 
and summarized in Fig. [7] by the time trend of the macroscopic and microscopic 
fluxes across the bottleneck. The model reproduces now the oscillations of the 
passing direction at the bottleneck described e.g., in [T2J [TJl [TS] . In more detail, 
starting from the initial condition depicted in Fig. [6^,, pedestrians of population 
2 are induced to stop at the bottleneck while those of population 1 go through 
at one side (Fig. ^jp, Fig. [7] for 4.5 < t < 8.5). After some time, population 2 
reorganizes and stops the flow of population 1 (Fig. |6j;, Fig. [7]for 8.5 < t < 11), 
then its larger mass stuck at the bottleneck gives it locally the necessary strength 
for repelling opposite walkers and gaining room in the middle (Fig. [6ji, Fig. 
for 11 < t < 15). Some walkers of population 1 remain trapped by the stream 
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Figure 8: Test 3. Alternate lanes through the bottleneck (left) and correspond- 
ing macroscopic fluxes (right) emerging with the fully microscopic dynamics. 
Notice that both fluxes are identically zero for t > 8 because by then both 
populations have completely flowed across the door 



of population 2 and cannot access the passage (Fig. [6^, Fig. [7] for t « 15) until 
most of population 2 has flowed through (Fig. [6]:, Fig. [7] for 15 < t < 17.5). 

From the modeling point of view, the difference with the case 6 = is that 
the multiscale coupling shifts some macroscopic mass (30% in this example) onto 
the microscopic pedestrians, all other parameters and initial conditions being 
unchanged. The inhomogeneous distribution of this microscopic mass induces 
a break of symmetry between the interfacing populations, which finally leads 
to an alternate occupancy of the passage according to the repulsion prevailing 
locally in space and time. 

Setting 6 = 1, which amounts to shifting the whole mass onto the micro- 
scopic pedestrians and having the microscopic scale lead the dynamics, produces 
instead the outcome displayed in Fig. [8j Now the microscopic granularity fully 
dominates, hence the stream is the most fluent one (cf. also Test 2). As a result, 
the bottleneck interferes less with the stream than in the previous cases, and 
the model reproduces the alternate oppositely walking lanes (Fig. [8^,) exten- 
sively observed as one of the main effects of self-organization in real crowds (cf. 
rr51 ITS] and references therein) . The time trend of the macroscopic flux across 
the bottleneck (Fig. ^p) confirms that the two populations flow simultaneously 
through the passage, with comparable fluxes, in the interval 1.8 < t < 8. After 
the time t = 8 the macroscopic fluxes are identically zero because by then both 
populations have completely flowed across the door. 

Test 4: Macroscopic effect of a microscopic leader 

In this last test we outline a capability of the multiscale model, not yet high- 
lighted so far, which will surely deserve further investigation. More precisely, we 
are referring to the use of the microscopic scale for modeling some features of the 
system which could not be described in a purely macroscopic framework, but 
which nonetheless affect the macroscopic dynamics. This is essentially different 
from the previous tests, where the same effects, such as repulsion and obstacle 
avoidance, were described at both scales. The main novelty here is that the 
system includes a microscopic term with no macroscopic counterpart. 

We consider the case of a crowd following a leader, for instance a group of 



MULTISCALE MODELING OF GRANULAR FLOWS 



26 



(a) 



(c) 



w 



0.4 
0.3 
0.2 
0.1 






1 - 



(b) 



(d) 



P 



0.4 
0.3 
0.2 
0.1 



0.4 
0.3 
0.2 
0.1 





Figure 9: Test 4. (a) initial condition, (b) the group assumes an elongated 
configuration, (c) the group is formed and follows the leader, (d) the leader 
waits for the group while the group moves on 



tourists and their guide. The leader is a microscopic pedestrian who behaves in 
a different way with respect to all of the other group members: s/he is the only 
one informed of the way to go, hence s/he walks with a pre-assigned velocity 
(0.4i in this example) independently of the others (i.e., s/he does not interact 
with the rest of the group). S/he only stops when her/his distance from the 
group becomes too large. The followers have zero desired velocity, because they 
are not informed of the way to go, and experience both frontal attraction and 
frontal repulsion with their group mates, including the leader. (Like in Test 1, 
the angle a xy in Eq. ( 20 1 is computed by assuming conventionally v^cs = i) • 
Attraction acts against group dispersion, and is needed especially in order for the 
crowd to follow the leader. Instead, repulsion is intended for collision avoidance 
among group mates. The radii R a , R r are equal (cf. Table [T]), in particular R a 
is so small that the tail of group does not feel the leader ahead. 

The group starts from the square-shaped distribution depicted in Fig. [9^,, 
with the leader in front. Then, after a transient (Fig. |9Jd), it assumes a hori- 
zontally elongated shape (Fig. [9J;) as a result of joint attractive and repulsive 
effects. With no leader such a configuration would be an equilibrium, as at- 
traction and repulsion balance. However, as soon as the leader starts moving 
forward undisturbed, pedestrians at the front, who can feel him directly, are 
attracted and move forward in turn. At the same time, pedestrians at the rear 
are attracted toward group mates in front. This makes the information on the 
way to go travel backward across the group, which ultimately moves forward as 
a whole. 

It is worth stressing again that at the macroscopic scale there is no counter- 
part of the microscopic leader. This implies that the macroscopic interaction 
velocity v[M t ] is not affected by the microscopic leader, therefore the macro- 
scopic mass feels the latter only through the microscopic interaction velocity 
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v[m t \. 

This test shows that our multiscale framework is suitable to reproduce a well 
known feature of self-organizing groups, namely the fact that a small number 
of informed agents can move the whole group in the desired direction [8]. In 
particular, thanks to the multiscale coupling, this effect is appreciable also at 
the macroscopic scale, which would not be suitable by itself to model differences 
among the individuals. 

7 Conclusions and future research 

In this paper we have presented a measure-based multiscale method for modeling 
pedestrian flow. We point out that neither a purely microscopic ODE-based 
approach nor a purely macroscopic PDE-based approach is new in the literature 
for this kind of application. The novelty here is the way of coupling the two scales 
in a rigorous mathematical framework. This is possible thanks to the measure- 
theoretic approach, which makes no a priori distinction between the scales, 
and to the fact that, by a proper scaling, the microscopic and the macroscopic 
model can reproduce the flow of the same mass of pedestrians with comparable 
outcomes (cf. the numerical test 1). We stress that introducing microscopic 
heterogeneity in a macroscopic model is not straightforward. Adding random 
disturbances to the macroscopic variables may lead to apparently good results 
but it cannot be mathematically nor physically justified in an averaged context. 
Instead, our method allows to add granularity to the macroscopic flow and to 
preserve at the same time physical meaning and mathematical rigor. 

From the modeling side, it is worth noticing that a macroscopic model of 
pedestrian flow is useful to get overall distributed information, especially in 
connection with design, control, and optimization issues. However, as demon- 
strated by our numerical simulations, a certain amount of granularity is often 
crucial to catch some aspects of self-organization in crowds triggered by the 
microscopic inhomogeneities of the flow (cf . the numerical test 3) . Of course, it 
has to be expected that the outcome at whatever scale partly depends on the 
tuning of the parameters of the model. Therefore, it is not our purpose to state 
that the multiscale approach is always better (i.e., more realistic) than either 
the microscopic or the macroscopic approach by itself. Rather we believe that 
the proposed technique offers a convenient way to make the two scales interact 
and jointly contribute to the final result. 

The present form of our multiscale approach is mainly concerned with the 
same mass of pedestrians modeled at both the microscopic and the macroscopic 
scale. The multiscale coupling is then realized by scale interpolation. However, 
the numerical test 4 demonstrates that the framework is suitable also to model 
features at either scale, which have no explicit counterpart at the other scale and 
nonetheless affect crucially the overall dynamics. As a research development, 
we plan to further generalize our multiscale approach in this direction, having 
in mind specific applications related to traffic flow. 

Pedestrians and cars share some relevant features, such as a desired velocity 
driving them toward specific destinations and a frontally restricted visual field. 
Actually car movements are much more constrained than pedestrians', hence 
self-organization is more limited, however not completely inhibited. For exam- 
ple, an application quite considered in the technical literature [HI [33] concerns 
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mixed traffic conditions with few mopcds within a flow of cars. 
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